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We present the temperature and layer dependent electronic structure of a 20-layer EuS(100)-filni 
using a combination of first-principles and model calculation, the latter based on the ferromag- 
netic Kondo-lattice. The calculated thickness-dependent Curie temperature agrees very well with 
experimental data. The projected 5d-bandstructure is at finite temperatures strongly influenced 
by damping effects due to spin exchange processes. Spin-split unoccupied 5d-surfaces states are 
found with a Stoner-like collapsing for increasing temperature towards the Curie point and with an 
exponential decay of spectral weight with increasing distance from the surface. 
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Europiumsulphide is a prototype ferromagnetic semi- 
conductor, crystallizing in the rocksalt structure with a 
lattice constant a = 5.95 A. The Eu'^^ ions occupy lattice 
sites of an fee structure. The just half-filled 4f shell of the 
rare earth ion creates a strictly localized magnetic mo- 
ment of li-j, B according to the ground state configuration 
^ 5*7/2 ■ The moments are exchange coupled resulting in 
a ferromagnetic order for temperatures below the Curie 
temperature Tc = 16.57K As to the purely mag- 
netic properties bulk-EuS is considered rather well under- 
stood being an almost ideal realization of the Heisenberg 
model, where the exchange integrals can be restricted to 
nearest [Ji/ks — 0.221K) and next nearest neighbors 
(J2AB = -O.IOOK), only 0. 

Recently, there has been an increasing interest in thin 
EuS films, in particular what concerns the thickness de- 
pendence of the ferromagnetic transition temperature jj| 
which is frequently discussed in terms of finite-size scaling 
(dimensionality crossover) Q . Investigations of this kind 
have already earlier been done for films of the metallic 
counterpart Gadolinium j^. In addition, Gd has pro- 
voked numerous research activities with respect to its 
extraordinary magnetic surface properties as, e. g., a 
possibly enhanced Curie temperature of the (0001) sur- 
face compared to that of bulk Gd . 

While the Heisenberg model provides an excellent de- 
scription of the purely magnetic properties of local- 
moment insulators such as EuS, EuO or metals like Gd it 
is of course inadequate for describing electronic and mag- 
netooptic effects. So it cannot explain the striking tem- 
perature dependence of the 5d conduction bands, first 
observed for the ferromagnetic europium compounds as 
red shift of the optical absorption edge for electronic 4f- 
5d transitions upon cooling below Tq 0. The reason is 
an interband exchange coupling of the excited 5d elec- 
tron to the localized 4f moments that transfers the tem- 
perature dependence of the ferromagnetic moment state 
to the unoccupied conduction band states. Recently the 
same temperature effects have been verified for film struc- 
tures in the case of EuO by spin-resolved x-ray absorption 



spectroscow I3 and by quasiparticle band structure cal- 
culations I3, ISj ■ In ref. [9j the theoretical prediction of 
a magnetic surface state has led to the speculation that 
the temperature dependence below Tq may give rise to 
a surface insulator-halfmetal transition accompanied by 
a huge magnetoresistance effect. It is one of the goals 
of our present study to teach whether this spectacular 
proposal may hold for EuS too. 

We present in this paper a study of the temperature 
and layer dependent electronic structure of thin EuS- 
films by combining a many-body model with a first- 
principles bandstructure calculation. The final goal are 
statements about magnetic stability in terms of the Curie 
temperature and the existence of surface states. 

In the following we assume EuS film structures with 
two surfaces parallel to the fcc(lOO) crystal plane and 
consisting altogether of d equivalent parallel layers. The 
lattice sites R^q. within the film are indicated by letters 
a, (3, ... denoting the layer, and by letters i, j, ... num- 
bering the sites within a given layer, for which we assume 
translational symmetry in two dimensions. Each lattice 
site is occupied by a magnetic moment, represented by 
a spin operator S^a , due to the half-filled 4f shell of the 
Eu^"*" ion. The exchange coupled moments (spins) are 
certainly well described by a Heisenberg Hamiltonian: 

As mentioned the exchange integrals J^j^ can be re- 
stricted to nearest ( Ji) and next nearest ( J2) neighbors. 
The dipolar energy of EuS films is taken into account by 
a single-ion anisotropy Dq. This helps to overcome the 
Mermin- Wagner theorem ^3 which forbids a collective 
moment order in films of finite thickness with isotropic 
Heisenberg exchange . Since the exchange constants are 
derived from a low-temperature spin wave analysis 0, 
we fix Dq, somewhat arbitrarily, by the requirement that 
our theory yields the experimental Tc for bulk EuS. We 
find Do/kB = 0.375K 

For the 5d conduction electrons we use the partial op- 
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erator H^d- 



^5d — Z^ijaf} l^7nm' -^ijaf3 ^iama'^j Pni' <J 



(2) 



Hama- {ciama) IS thc Creation (annihilation) operator of 



an electron with spin a at site Rm in the orbital m. The 
Tj™^ describe the electron hopping from Rj^ to Rm 
with a possible orbital change (m' — » m). We require 
that these single-electron energies do not only account 
for the kinetic energy and the influence of the lattice 
potential, but also for all those interactions which are not 
explicitly covered by the model Hamiltonian. Therefore 
we take them from an LDA calculation. Details are given 
below. 

Many important properties of local-moment materials 
such as EuS can be traced back to thc interaction between 
the 4f and 5d partial systems. Starting from the very gen- 
eral on-site Coulomb interaction between electrons of dif- 
ferent orbitals it can be shown |^ 0| for the special case 
of EuS (half-filled spin-polarized 4f shell, empty conduc- 
tion bands) that the interband interaction can be written 
as: 



(3) 



J is the corresponding exchange coupling constant. Fur- 
thermore, we have used the standard abbreviations: 

riimaa = Araaa^^raaa; ^ S"^ ± iS^ . The firSt term 

in (j2Jl describes an Ising-like interaction between the lo- 
calized 4f spin and the spin of the itinerant 5d electron, 
while the two others provide spin exchange processes be- 
tween the two subsystems. Spin exchange may happen 
in three different elementary processes: magnon emis- 
sion by an itinerant J, electron, magnon absorption by 
an I electron and formation of a quasiparticle {magnetic 
polaron). The latter can be understood as a propagating 
electron dressed by a virtual cloud of repeatedly emitted 
and reabsorbed magnons corresponding to a polarization 
of the localized-spin neighborhood. 

We believe that the total model Hamiltonian 



H = K 



4f 



5d 



H, 



4f-5d 



(4) 



incorporates the main physics of the ferromagnetic local- 
moment insulator EuS. It can be considered as the multi- 
band version of the ferromagnetic Kondo-lattice model 
(KLM). To get a realistic description of EuS we try to 
cover all those interactions, which do not explicitly ap- 
pear in our model Hamiltonian, by a proper renormal- 
ization of the single-particle energies. For this reason we 
performed a band-structures calculation using the tight- 
binding LMTO-atomic sphere approximation program of 
Andersen 0, ■ Difficulties typical of LDA arise with 
the localized character of the 4f levels. A "normal" LDA 
calculation produces wrong 4f positions. To circumvent 
the problem we consider all the seven 4f electrons as core 



electrons in one spin channel, since the 4f levels enter our 
study only as localized spins in the sense of H4{ in 
We have to choose the proper single-particle input in such 
a way, that a double counting of just the decisive inter- 
band exchange explicitly by the model Hamiltonian 
Q and implicitly by the LMTO input, is avoided. The 
most direct solution of this problem would be to switch 
off the interband exchange H^f.^^ in the LDA code, what 
turns out to be impossible. We can exploit the fact that 
the non-trivial many-body problem of the KLM is ex- 
actly solvable for a ferromagnetically saturated semicon- 
ductor, e. g. EuS at T = OK The result is especially 
simple for the t spectrum because the "f electron cannot 
exchange its spin with the fully aligned 4f spins. Only 
the Ising term in iJ4f_5d takes care for a rigid shift of 
the total t spectrum by ^JS. Thus we can use without 
any manipulation the | dispersions of the LMTO calcu- 
lation, which by definition works for T = OK, for the 
EuS-hopping integrals in There is no need to switch 
off i/4f-5d because in this special case it leads only to an 
unimportant rigid shift. It should be stressed, however, 
that the | spectrum, on the contrary, is strongly influ- 
enced by the interband exchange, even at T = OK where 
it is still exactly calculable [13 (see fig. 4). 

Because of the absence of conduction electrons in the 
semiconductor EuS, at least at low temperatures, the 
magnetic ordering of the localized 4f moments will be 
unaffected by the band states, because an RKKY-type 
contribution cannot occur. On the other hand the rel- 
evant superexchange mechanism is hidden in thc ex- 
change integrals J°'j^ in Q), which are taken from the 
experiment (2| The magnetic part of the problem can 
therefore be solved separately via the extended Heisen- 
berg Hamiltonian Hif. We employ a generalized version 
of the well-known Tyablikow approach, from which we 
know that it yields convincing results in the low- as well 
as high-temperature region. Details of the method are 
presented in 0. The result is a distinctly layer- and 
temperature-dependent magnetization of the 4f moment 
system. The Curie temperature of the EuS-film exhibits 
a strong thickness-dependence, starting at about 2K for 
the monolayer and approaching the bulk-value for more 
than 25 monolayers (Fig. The Tq of the 20 layer 
film, we are going to discuss below, amounts to 16.28K. 
We find a remarkable agreement of our theoretical re- 
sults with experimental data found by Stachow-Wojcik 
et al. 0. A similar Tq behavior, derived from suscep- 
tibility measurements, is reported by Farle et al. for 
the metallic counterpart Gd. The decrease of Tq with 
decreasing film thickness is understood as typical finite 
size scaling accompanied by a dimensionality cross-over. 

To obtain the temperature-dependent electronic struc- 
ture of a EuS(lOO) film the T = O-f dispersions from 
our tight binding-LMTO band calculation are needed as 
single-particle input for the model evaluation. To ac- 
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FIG. 1: Curie temperature as function of film thickness (num- 
ber of EuS monolayers). Our theoretical results: o - - o. The 



experimental data are taken from ref. 



16.57K. 



count for the film geometry, we define a supercell con- 
sisting of n consecutive EuS(lOO) layers followed by m 
layers of empty spheres, i. e. periodically stacked EuS n- 
layer films, isolated from each other by m layers of empty 
spheres. It turned out that m = 5 is large enough to 
guarantee truly isolated EuS films. We do not consider 
surface relaxation and reconstruction, which might be 
present in EuS films. 

Similarly to the magnetic part the electronic part 
can be treated separately because magnon energies are 
smaller by some orders of magnitude compared to other 
energy terms as the exchange coupling constant J or the 
conduction band width. So we can disregard the mo- 
ment operator H^f when calculating the electronic self 
energy. That does not at all mean that the magnetic 
4f moments do not influence the electronic quasiparticle 
spectrum. That is rather done by local-spin correlations 
such as {S'), (S^S^), {{S')^),- ■ •, which are to a great 
extent responsible for the temperature-dependence of the 
electronic spectrum and have to be derived from H^f. 
To get the electronic self energy we apply the moment- 
conserving decoupling approximation (MCDA) for suit- 
ably defined Green functions 0, . In the set of equa- 
tions which constitute the formal solution for the elec- 
tronic self energy there appear the just mentioned 4f- 
spin correlations. Although the MCDA is partly based 
on an intuitive ansatz, it has been proved to be quite a 
reliable approach to the sophisticated many-body prob- 
lem of the multiband-KLM. Interpolating self energy ap- 
proaches [l^ . which fulfill a maximum number of exact 
limiting cases, as well as a systematic projection operator 
method yield essentially identical results. We do not 
resent details of the MCDA but refer the reader to ref. 

In summary our theory contains only one parameter, 
namely the exchange coupling J. We assume that it is 
the same value for the film as for bulk-EuS, J = 0.23eV 
[m . Note that J is not really a free parameter but being 
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FIG. 2: Projected unoccupied (5d) t band structure (T = OK) 
for a 20 layer EuS (100) film. The shaded region belongs to the 
respective projected band structure of bulk EuS. The vertical 
broken line identifies the position in k-space of the surface 
state near the energy zero, which is investigated in Fig. |3 
The energy zero has been chosen arbitrarily. 

Fig. 121 shows the unoccupied (5d) t band structure of 
a 20 layer EuS(lOO) film at T = OK, calculated with 
the tight-binding-LMTO method. These results repre- 
sent the single-particle input for our model calculation. 
Note that the energy zero has been chosen arbitrarily, not 
at all coinciding with the Fermi edge. For simplicity we 
restrict our presentation to the 5d states only For 
comparison we also indicate in Fig. |2] the corresponding 
projected (5d) bulk-band structure, which has the den- 
sity of states shown in the upper part of Fig. 2 in 
Above all that helps clearly to identify surface states. A 
surface state is a state which appears in the forbidden 
region where no bulk states occur. Another character- 
istic feature of a surface state is the exponential decay 
of its spectral weight with increasing distance from the 
surface. To demonstrate this, we have plotted in Fig. |3| 
the up-spin spectral density as function of energy for the 
cut, which is marked in Fig. |21 by a broken line, i. e. 
for the wave vector at |rX. Furthermore, the spectral 
density is represented for different layers of the 20-layer 
film, a = 1 denotes the surface layer and a = 10 the 
middle layer. Note that the spectral density consists for 
T = and cr =| exclusively of delta- functions represent- 
ing quasiparticles of infinite lifetimes. Only to visualize 
the spectrum we have added a small imaginary part to 
the self energy getting therewith spectral density peaks 
of finite widths. The structure within the shaded region 
is obviously a surface state. Its spectral weight decreases 
exponentially with increasing distance of the respective 
layer from the surface. Another surface state appears 
at the bottom of the spectrum. We interpret this surface 
state as the analogue to that of EuO which has been spec- 
ulated in ref. to give rise to a surface metal-insulator 
transition when cooling down from Tc to T = 0. Fig.^ 
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FIG. 3: Spectral density {T — 0, a =|) as function of energy 
for the wave vector at §rX (dashed Une in Fig. |5J and for 
different layers {a = 1,2,3, 10) of a 20-layer EuS(lOO) film. 
The inset shows on a logarithmic scale the spectral weight of 
the state within the grey column for the surface and the two 
layers next to the surface, a — 1(10) means surface (middle) 
layer. 



shows the temperature dependence of the spectral den- 
sity, again for the wave vector at |rX to include the 
above-discussed surface state. At T = (ferromagnetic 
saturation) the spiky structure of the t spectral density 
refers to quasiparticles of infinite Hfetimes. The t elec- 
tron has no chance to exchange its spin with the par- 
allel aligned localized 4f spins. It therefore propagates 
through the lattice without any scattering. On the other 
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FIG. 4: Spectral density of a 20-layer EuS(lOQ) film at frX 
as function of energy and that for three different temperatures 
(T = 0, O.SSTc, Tc). Upper halves for the up-spin spectrum, 
lower halves for the down-spin spectrum. Full lines: surface 
layer (a = 1), thin lines (shaded regions): middle laver {a — 
10). Tc = 16.28K for the 20-layer film. The vertical broken 
lines indicate the position of the spin-split surface state; left 
line for cr =|, right line for cr =J,. 



hand, a J, electron can be scattered by magnon emission 
or can form a magnetic polaron what gives rise to a dra- 
matic lifetime broadening of the spectrum even at T = 0. 
We want to stress that the (T = 0, cr =J,) calculation is 
exact: For finite temperature, i. e. finite magnon den- 
sities, the t electron, too, can be scattered by magnon 
absorption with a concomitant spin flip. Such processes 
lead to a rather structureless 5d excitation spectrum, 
which, nevertheless, carries a distinct layer dependence. 
The surface state, identified in Fig. 13 is spin split at 
T = OK by some 0.8eV due to the exchange coupling to 
the ferromagnetic 4f moment system. The induced ex- 
change splitting is strongly temperature-dependent with 
a Stoner-like collapsing for T Tc- The same behavior 
is shown up by the surface state at the bottom of the 
spectrum. An analogous feature has been reported for 
EuO in ref. . This distinct temperature dependence in 
the ferromagnetic phase appears at first glance somewhat 
astonishing because it happens to the unoccupied and "a 
priori" uncorrelated 5d energy states of semiconducting 
EuS. 

By a combination of an LDA-band structure calcu- 
lation with a many-body evaluation of the multiband 
Kondo-lattice model a pronounced temperature depen- 
dence of the (empty) 5d conduction bands of EuS(lOO) 
films could be demonstrated. Exchange split magnetic 
surface states show up a Stoner-like collapsing for T — > 
Tc- The Curie temperature exhibits a characteristic film- 
thickness dependence understandable as a typical finite 
size effect. 
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